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Abstract 

By analytically solving some simple models of phase-ordering kinetics, 
we suggest a mechanism for the onset of non-equilibrium behaviour in 
colloid-polymer mixtures. These mixtures can function as models of 
atomic systems; their physics therefore impinges on many areas of ther- 
modynamics and phase-ordering. An exact solution is found for the mo- 
tion of a single, planar interface separating a growing phase of uniform 
high density from a supersaturated low density phase, whose diffusive 
depletion drives the interfacial motion. In addition, an approximate so- 
lution is found for the one- dimensional evolution of two interfaces, sepa- 
rated by a slab of a metastable phase at intermediate density. The theory 
predicts a critical supersaturation of the low-density phase, above which 
the two interfaces become unbound and the metastable phase grows ad 
infinitum. The growth of the stable phase is suppressed in this regime. 

PACS numbers: 82.70.Dd, 64.60.My, 05.70.Ln 

1 Introduction 

Metastability and phase transition kinetics in condensed matter are rapidly ex- 
panding fields of research [|T], HI. Experimentally, it is increasingly realised that 
colloidal systems exhibit unique properties which make them attractive candidates 
for studying these topics. Advances in synthetic chemistry mean that suspensions 
of particles of well-characterized shapes and sizes, and with precisely tunable inter- 
particle interactions, can be routinely produced. The equilibrium phase behaviour 
of such model colloids can be studied experimentally and understood in some detail 
using the standard tools of statistical mechanics |§. Moreover, colloids are much 
larger than atoms, the typical dimension of a colloid being L pa 0.5/im, and the 
solids they form (crystal, glass, gel) are mechanically weak, a typical modulus being 
G ~ ksT/L 3 pa 10~ 2 Nm~ 2 . Thus, a colloidal crystal can be 'shear melted' to the 
metastable colloidal fluid state simply by shaking. Structural dynamics of colloidal 
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systems are also slow, the scale being set by the time a free particle takes to dif- 
fuse its own dimension, r# ~ L 2 /D. Estimating the diffusion coefficient D by the 
Stokes-Einstein relation for a spherical particle of radius L ~ 0.5/zm suspended in 
a liquid of viscosity rj ~ 10~ 3 Nm~ 2 s, we get D = k^T '/QnrjL ~4x I0~ 13 m 2 s~ 1 , so 
that the characteristic relaxation time in colloidal systems is tr > 1 s. Thus, the 
crystallization of a metastable colloidal fluid can take minutes, hours or even days. 
These two features, the ease with which metastable states can be prepared and the 
long relaxation times involved, mean that the study of phase transition kinetics and 
metastability in model colloids is gaining increasing attention. For example, there 
is a growing literature on crystallization kinetics in nearly-hard-sphere PMMA col- 
loids [|J. The glass transition in the same system is now recognized as being one 
of the simplest examples of that phenomenon, and has become a testing ground for 
sophisticated theories |J. 

In recent work aimed at elucidating the equilibrium phase behaviour of mixtures 
of nearly-hard-sphere colloids and non-adsorbing polymers, we have also observed a 
rich variety of non-equilibrium behaviour |], [7| . The well-characterized nature of our 
experimental system, coupled with detailed data from small-angle light scattering 
[0], have allowed us to suggest a possible link between the onset of formation of non- 
equilibrium phases and the presence of a 'hidden', metastable minimum in the free 
energy landscape. In this paper, we review the experimental evidence leading to this 
suggestion, and model the kinetic consequences of such a three-minima free energy 
landscape in one dimension using a continuum approach. We conclude by pointing 
to the possible applications of this model in other 'soft' systems. Some of our work 
relies on results from a companion paper ||; the main ideas are summarized in Q. 

2 Non-Equilibrium Behaviour in Colloid- Polymer 
Mixtures 

Free polymer alters the phase behaviour of colloids via the depletion mechanism ]I(J . 
Exclusion of polymer from the region between the surfaces of two nearby colloidal 
particles gives rise to an unbalanced osmotic pressure pushing the particles together, 
resulting in an attractive depletion potential, U^ ep . In the case of hard-sphere-like 
colloids, the total interparticle interaction [|ll|] in the presence of polymers is given 
approximately by 




for r < a 

for a < r < o + 2r g 
for r > a + 2r g 
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where a = 2a is the particle diameter and U p (fj, p ) is the osmotic pressure of the 
polymer as a function of the polymer chemical potential, /i p . Kveriap is the volume 
of the overlapping depletion zones between two particles at an inter-centre separation 
of r. Explicitly 
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where £ = r g /a denotes the ratio of the radius of gyration of the polymer, r g , to the 
radius of the colloidal particle, a. 

Theory |10| predicts and experiments show that the phase diagram of a mix- 
ture of hard-sphere colloids and non-adsorbing polymer is a function of the size ratio, 
£. The phase diagram at large enough £ is reminiscent of that of a simple atomic 
substance, figure |l|a, with the polymer chemical potential (/x p ) playing the role of 
an inverse temperature. (Temperature itself is not an axis of the phase diagram 
because, in systems with only excluded-volume interactions, no energy scale exists 
O.) Colloidal gas, liquid and crystal phases are possible, with regions of binary 
coexistence between pairs of these phases and triple coexistence at a particular value 
of fjL p . As the size ratio £ is decreased, the region of colloidal gas-liquid coexistence 
shrinks until, below a critical value £ crit , the gas-liquid coexistence region disappears 
altogether from the equilibrium phase diagram, figure |l|b; now, only colloidal fluid 
or crystal phases exist in equilibrium. Experimentally, £ cr i t ~ 0.25. 

Focus now on the phase diagram for the case of £ < £ cr i t , figure Experiments 
confirm that moderate concentrations of a small non-adsorbing polymer cause a sus- 
pension of hard spheres to phase separate into coexisting colloidal fluid and crystal 
phases. At higher colloidal concentrations, or higher fi p , across a non- equilibrium 
boundary (NEB), crystallization is suppressed 0, 0. Instead, a variety of non- 
equilibrium aggregation behaviour is observed. Detailed small-angle light scattering 
studies have allowed the classification of the kinetic behaviour above the NEB 
into three types. Just across the NEB, the behaviour is 'nucleation like'. This is 
characterized by an initial latency period consistent with nucleation, after which 
dense, amorphous droplets separate out, forming an amorphous sediment which be- 
gins to crystallize only on a much longer time-scale. The 'nucleation-like' regime is 
of primary interest to us in this article. At higher or higher colloid concentra- 
tion, the behaviour becomes 'spinodal-like', switching finally to a 'transient gelation' 
regime for the densest systems. 

It should be stressed that the experimentally observed NEB is sharp and highly 
reproducible. In fact, it has been suggested by one of us Q that this boundary is 
given by a 'hidden' gas-liquid binodal. Following the theoretical approach in JlQ| . 
we write the free energy, F(p, fi p ), of a colloid-polymer mixture as a function of the 
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Figure 1: Phase diagrams of colloid-polymer mixtures in the colloid concentration 
(p)-polymer chemical potential (p p ) plane, for size ratio £ (a) above and (b) below 
£ cr i t . The dotted curve is a "hidden" gas-liquid binodal. For (b), the free energy 
plots for two polymer chemical potentials, /zp° w and yUp" gh , are shown in figure ^]a and 
b respectively. The values of p for the various boundaries at p = p^ gh are given by 
the common tangent construction in figure 0b. 



colloid volume fraction (or 'density'), p, and the polymer chemical potential, p p . F 
can be calculated within a mean-field framework for a disordered arrangement of 
colloids and polymers, the 'fluid branch', and an ordered arrangement of colloids 
with polymers randomly dispersed, the 'crystal branch'. At low polymer chemical 
potentials, the fluid and crystal branches each show a single minimum (figure ||]a). 
This gives rise to single-phase fluid (p < p g ), fluid-crystal coexistence (p g < p < p s ), 
or single-phase crystal (p > p s ). The colloid concentrations in coexisting fluid and 
crystal phases are obtained by the common tangent construction (see, e.g., |13| ). 

At higher polymer chemical potentials, however, the fluid branch of the free en- 
ergy shows a 'double minimum' structure (figure |2|b). At larger polymer to colloid 
size ratios, this double minimum can give rise to a region of colloidal gas-liquid co- 
existence in the equilibrium phase diagram, figure |I]a |10[ [11]. For small polymers, 
however, the theory predicts only separation into fluid and crystal phases. Never- 
theless, the 'metastable gas-liquid binodal', which is 'buried' within the fluid-crystal 
coexistence region predicted by equilibrium thermodynamics, can still be traced out 
(figure |l]b). The equilibrium phase boundary and the metastable gas-liquid binodal, 
calculated using the theory in |T0| , compare well with the experimental equilibrium 
fluid-crystal coexistence boundary and the non-equilibrium boundary respectively 
[0], giving rise to the speculation that the suppression of crystallization, which is 
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(a) (b) 




Figure 2: Fluid and crystal branches of the free energy density versus colloid volume 
fraction (a) at low polymer chemical potential (e.g. /4 ow in figure |]b), when the fluid 
branch has a single minimum, and (b) at higher polymer chemical potential (e.g. 
fip lgh in figure |l]b), when two minima exist in the fluid branch, corresponding to gas 
and metastable liquid phases. The common tangent construction is used to find 
the densities of the coexisting equilibrium gas and solid phases (p g ,p s ), and the 
metastable gas-liquid and liquid-solid binodal points (pa,Pb), (pcPd)', see figure 
0b. 



the predicted, equilibrium phase behaviour, is linked to the presence of the hidden 
gas-liquid binodal. 

In this and a preceding, companion paper [§], we provide a theoretical basis for 
this speculation. We model the effects of a metastable third minimum in the free 
energy curve on the kinetics of phase ordering, focusing on post-nucleation kinetics. 
A system which has been homogenized (e.g. by shear-melting) to some uniform 
density between p g and p s (see figure §) must evolve towards the equilibrium state, 
by separating into regions of density p g and p s . If f(p) has positive curvature 
at the initial density, then the evolution begins by nucleation, whereby a thermal 
fluctuation causes part of the system to cross the barrier between the two stable 
wells in /. Once such a nucleus has formed (in the early stage of phase ordering), it 
grows by diffusively depleting the surrounding, supersaturated medium (the growth 
stage [fUJl). 

In the previous paper ||, we consider the implications for short range (on the 
length scale of the transition zone between neighbouring phases) density variations 
and stability, in the context of quasi-steady-state motion. There, and in Ref. ||, 
it is shown that an interface between regions of different densities may split, in the 
presence of a metastable minimum, into two parts. The purpose of the present 
article is to model the evolution of a split interface during the growth stage, consid- 
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ering large length scales (so interfaces appear as sharp steps, and relatively gentle 
concentration gradients control the diffusive dynamics) and time-scales (so motion 
is time-dependent, not steady-state). We will find that growth dynamics involving 
split interfaces contrasts greatly with normal growth behaviour. 

The plan of this paper is as follows. The model introduced in section |3], involv- 
ing a nucleus with an ordinary, unsplit interface, sets the scene for the methods 
employed. It is a completely soluble model for the motion of a single domain wall. 
In section |J we analyse a more elaborate scenario, where two domain walls (the two 
parts of a "split" interface) enclose a region of metastable density, and compete for 
the flux of condensing material. This model exhibits behaviour consistent with the 
non-equilibrium boundary in the colloid-polymer phase diagram. The limitations of 
our work are discussed in section |J while applications to related systems are given 
in section |6|. 

In the models presented, growth is limited by diffusion of matter only. This 
simplification holds in many soft condensed systems, for which diffusion of latent 
heat is irrelevant since the entropy is dominated by the degrees of freedom of the 
solvent EJ. Moreover, we assume that the system may be described by a single 



conserved order parameter (such as the density of colloidal particles). Formally this 
does not apply in colloid-polymer mixtures, where the polymer concentration is a 
second conserved quantity. However, it should be a reasonable model for small £, 
such that the polymers diffuse much more rapidly than the colloids. Non-conserved 
order parameters describing crystallinity are also omitted, under the assumption 
that such variables relax quickly compared with the diffusive order parameter, and 
are therefore not rate-limiting. 

3 Exact Solution of Diffusion Equation for Two- 
Phase Separation 

In this section, we find an exact solution of a problem pertinent to the intermediate- 
to- late-stage evolution of a two-phase system. It is well known [|14[] that, when a 
dense droplet grows into a surrounding, supersaturated medium, the speed of motion 
of the interface is proportional to t -1 / 2 . We solve an idealized one- dimensional model 
for the motion of a single interface in such a situation, and find the constant of 
proportionality as a function of the supersaturation of the ambient medium. The 
results obtained in this section will be used in section |] to investigate the effect on 
interfacial motion of a 'metastable minimum' in the free energy landscape. 
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3.1 Idealized Model 



Consider an infinite, one-dimensional system whose equilibrium state is two-phase 
coexistence, which has been homogenized to a uniform density (e.g. by shear melt- 
ing). In the system, a region of the high-density phase has nucleated at its equilib- 
rium density, and the interface between the high- and low-density phases has locally 
equilibrated. The domain wall proceeds to move by depleting the supersaturated 
low-density phase. Let the interface be located at x — x\(t) and let the density 
of the high-density phase be uniform at p — p s . So the evolution of the system is 
solely due to the dynamics of the low-density region. In this phase, the density is 
p g at the interface (x = Xi), and tends asymptotically to the supersaturated value 
p g + cr as x — > oo. The assumption that the density (and hence chemical potential) 
is fixed at the equilibrium value (given by the double tangent construction in figure 
0a) at the interface is physically valid if the interface moves sufficiently slowly ||. 
The problem of finding x\(t) is that of solving the diffusion equation 

dp = D d 2 p 
dt dx 2 

with boundary conditions p(oo,t) = p g + a and 

p(xi,t) = p g (3) 

with a moving boundary x\(t), whose velocity depends on the flux at that point 
since, by conservation of matter at the interface, 

(p s -Pg)x l =Dp'(x 1 ). (4) 

The initial conditions are xi = and p(x>x±) = p g + a, with a negative delta- 
function singularity at x = xi, since p(x±) = p g V t. The system is depicted, at 
some later time t, (and the initial condition is inset) in figure |3|. 

3.2 Solution 

Solving the diffusion equation with a moving boundary (a Stefan problem) is dif- 
ficult. The problem is overcome by replacing the semi-infinite region in which the 
diffusion equation is to be applied, by an infinite region, with a source or sink at 
Xi, whose strength s(t) is such that the region x > x\ cannot distinguish it from a 
moving interface. The behaviour in x < xi, which is just a construct of the method, 
should simply be ignored. The initial conditions in this region may be freely chosen 
to facilitate the solution. Let p (x^x±(0) , t — 0) = p g + cr and let y be defined as the 



7 



deviation from this ambient density; y = p — p g — a. So the equation to be solved 
is the diffusion equation with a source at x — Xi(t): 

% = D§-S( X -Mt))s(t) (5) 

with the wall position x\(t) and the construct s(t) fixed by the boundary conditions. 
Equation |5] is solved by 

y(x,t) = — dt'l dx' G{x — x' ,t — t') 8(x' — x\(t')) s(t') 

JO J-oo 



using the Green's function 



1 ( -x T 

G(x, t) = —j=^ exp 



Hence the density field is given in terms of s(t) and Xi(t) by 

/•* s(t') (x-xAt')) 2 
p(x, t) = p, + a- / dt' v ; = exp ; n '! . 



Applying equations |3] and |] gives us two integral equations in the two unknown 
functions Xi(t) and s(t), valid for all positive t; 



_ 

(xi(t) - Xl (t')) s(t) (x x (t) - x x {t')f _ n l / 2 (p s -p g ) , 

d (4D(t-t')) 3/2 P 4£>(t-tf) " 2D Xl[tl {0 

The right-hand side of equation || is independent of t, and hence the t dependence 
must be removed from the exponential in the integrand. This requires 



Xl (t) = aVDt, (8) 

and hence s(t) = by/(D/t), where a and b are constants in time. Substituting 
into equations || and ^ and evaluating the integrals (see appendix [A]) results in a 
closed-form expression relating the coefficient a to the relative supersaturation; 
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°— = ^ae^erfc§ (9) 



where erfc is the complementary error function, erfcx = 1 — erfx. The coeffi- 
cient b is not of particular interest, but we note that it is of the form (p s — p g ) x 
(function of — - — ). 
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Figure 3: Graph of density p against position x. On the left is the dense region 
which has been nucleated, and is growing at the uniform density p s , which is its 
equilibrium density. The interface is located at x — X\, which is a function of time, 
since matter is flowing down the concentration gradient in the sparse phase, and 
condensing onto the high-density region. At the base of the interface, the density 
is that of the coexisting gas, p g . The 'ambient' density far from the interface is 
p g + a. Inset: The initial (t = 0) configuration of the system, with a delta-function 
singularity at the base of the domain wall. 



We have chosen to use the initial condition £i(0) = in deriving equation but 
clearly, the origin of x\ is arbitrary, since the physics is translation- invariant, and 
does not depend on the initial size of the dense region. The velocity of the interface, 
on the other hand, is well defined: 

a ID . . 

*' = 2VT- (10) 

Equation [10] lends a physical meaning to the coefficient a, and so we will refer to it 
as the 'velocity coefficient'. 

Equation |9| is plotted in figure £|, for a as a function of a/ (p s — p g ). The validity 
of the model extends to negative supersaturations, which result in evaporation of 
the dense phase, and hence negative a. The curve plotted in figure [| has no special 
features in the negative quadrant. As the relative supersaturation tends to negative 
infinity, the behaviour of the velocity coefficient is given by 

a -> -2, In ( ~° ) . 

\ \Ps-pgJ 
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Figure 4: The velocity coefficient versus relative supersaturation. 



3.3 Linear and Non-Linear Regimes 

For small relative supersaturation (close to the origin of figure f|), equation || tends 
to a linear relation: 

2 f ^-\. (in 



"K \p s - Pg) 

From equation f|, we see that the flux onto the interface is j CO nd. = (p s ~ Pg) If we 



substitute equations pi] and 11 into this expression, we see that the flux of material 



from the dilute phase condensing onto the interface is 



jcond. « o-^JD/wt (12) 

in the linear regime. This is independent of {p s — p g ), and depends only on properties 
of the dilute phase. So the condensation flux becomes independent of properties of 
the interface in this regime. The physics behind this statement is as follows. The 
density at the base of the interface is maintained at the constant value p g and 
this gives rise to a gradient, and hence a flux, in the supersaturated dilute region. 
This gradient diminishes with time as material is depleted from the region, but is 
enhanced by the motion of the wall. It is clear from figure |3] that moving the wall to 
the right must accentuate the gradient. In the linear regime, the motion of the wall 
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is slow enough for this enhancement of the gradient to be insignificant, and the flux 
onto the interface varies with time as if the boundary to the diffusive region were 
fixed. In figure f|, we see that the model gives rise to a divergence (albeit unphysical) 
of the velocity coefficient at a relative supersaturation of unity. The reason for this 
is now clear. At this supersaturation, the non-linearity described above, whereby 
the interface motion enhances the gradient, becomes extreme. The interface can 
never deplete the 'low' density region, and hence the gradient (and therefore the 
flux) at its base remains infinite. 

Based on the exact solution of this simple model for interface motion (specifically 



equation 12), we proceed to study a more elaborate situation. 



4 Metastable Phase Evolution between Two Com- 
peting Interfaces 

It was noted in section ^ that the onset of non-equilibrium behaviour in colloid- 
polymer mixtures seems to be connected with the appearance of a third minimum 
in the free energy density, corresponding to a metastable liquid phase. Specifically, 
the formation of amorphous, non-equilibrium material from nucleation-like dynamics 
begins to occur at a supersaturation of the gas phase close to the hidden gas-liquid 
binodal. We now develop a model for the phase ordering dynamics in systems with 
a third minimum of this kind. 

In a system with such a three- well potential (like that illustrated in figure ^o), 
three different species of interfaces may exist between regions of the various locally 
stable densities, during intermediate- to late-stage ordering. If the gas, metastable 
liquid and solid phases are denoted g, I and s respectively, then g-l , l-s and g-s 
interfaces may move through the system. If we use the 'equilibrium interface' ap- 
proximation, as applied in section |3|, defining the densities above and below each 
interface to be the coexistence values for equilibrium between the two neighbouring 
phases, then a g-s interface has the fixed densities (p g ,p s ), given by the double- 
tangent construction [DJ| in figure 0b. Accordingly, at a g-l interface the densities 



are (pa,Pb), and at a l-s interface, (pc,Pd), which are the metastable coexisting 
values given in figure |2|b. 

Consider the idealized system introduced in section ||], where a region of the 
dense solid has appeared in an otherwise homogeneous, supersaturated system. The 
ambient density is within the gas well of the bulk free energy density. As before, 
the dense region will grow by depleting the supersaturated dilute region, but we can 
now imagine two alternative ways in which the growth can proceed. Either a g-s 
interface propagates into the dilute phase, and the system evolves as in section BL or 
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g-l and l-s interfaces form, and propagate separately; their motion being controlled 
by diffusion, both in the ambient dilute region and in the intervening region of 
metastable liquid. Once formed, a g-s interface propagates stably with respect to 
small perturbations || and cannot easily be split into g-l and l-s parts. Hence, which 
of these two modes of growth the system exhibits depends on which was initiated at 
the nucleation stage. The criteria for which mode is initiated during nucleation in a 
given system are unknown at present, although there has been some conjecture ||, |§. 
Let us accept that the split-interface (g-l and l-s) mode of evolution has begun in 
our one-dimensional system, and calculate the subsequent growth dynamics. 



4.1 Model System 

The density profile of the system in question, at some time t, is depicted in figure 
||. The g-l interface is positioned at xi(t), and the constant densities pA and ps 
immediately adjacent to it are marked, as are the fluxes ja{P) and Js(^) in and out 
of the interface, which are defined in the direction of the arrows, and are functions 
of time. The l-s interface is at x^(t). The flux of material of density pc into 
this interface is jc(t). There is no output flux, as the solid phase has a uniform 
density po- Let us make the model as general and physically realistic as possible 
by allowing different diffusion constants in the different phases: Di in the gaseous 
phase, and D 2 in the metastable liquid. The supersaturated density at infinity is 
defined as pa + &' in the figure. The 'adjusted supersaturation' a' is distinct from 
the supersaturation a used in section |3|. If the ambient (asymptotic) density is 
then the supersaturation is the deviation of this density from the equilibrium value 
(the stable binodal), which is p g in this system, i.e. a = p^ — p g . On the other 
hand, the adjusted supersaturation is the deviation from the metastable binodal, 
c' = Poo - Pa- 

Clearly, conservation of matter at the interfaces leads to the equations 

(Pb-Pa)xi = 3a- 3b (13) 
and (p D - pc)x 2 = 3c- (14) 



Using these relations, we could proceed in the same manner as in section |3]2|, but 
this time solving two coupled diffusion equations. The diffusion equation for the 
gaseous region would have one source, and that for the liquid region would have 
two. The strengths and positions of the sources would be fixed by equations IT3 and 



14] , and by fixing the constant densities pa, Pb and pc- Equation [L3| would couple 
the two equations. Proceeding in this manner, to try to find an exact solution for 
the evolution of the system would lead to five coupled integro-differential equations 
in five unknown functions (x±(t), x 2 (t) and the strengths of three sources). It is 
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Figure 5: Graph of density p against position x. The interface at x\ separates 
the gas phase from the metastable liquid, and is locally equilibrated so that the 
discontinuity in the density is between the metastable binodal values (pa,Pb)- The 
flux into (out of) the interface is ja Ub)- The flux j c of material of density pc 
impinges on the interface at x%, which separates the metastable liquid from the solid 
region of uniform density pr>. The 'ambient' {i.e. asymptotic) density of the gaseous 
phase is the supersaturated value pa + er'. Inset: The initial configuration of the 
system, with a uniform gradient between the domain walls, and a delta-function 
singularity in the gas phase at x%. 



not easy to spot the solutions to this system of equations, as it was for equations 
|6] and [7[ Instead we will make progress by introducing some physically reasonable 
appr oximat ions . 



4.2 Approximations 

The first approximation is to decouple the g-l interface from the gaseous driving 
region by using equation [12] for the input current: 



JA w (/y/DJrt. (15) 

This becomes exact as a' — > 0, and will presumably give qualitatively meaningful 
results at higher supersaturations, although with a' becoming some effective super- 
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saturation, deviating from the true value. We use a 1 rather than a in equation [0| 
as the boundary condition p{x\) — p A = was used in its derivation. 
As a second approximation, let 

r-. PB — PC , 

3b~3c~D 2 , 16 

X\ — %2 

which says that the gradient in the liquid region is approximately uniform. Intu- 
itively, this seems to be a reasonable assumption. If we imagine that, at some time 
into the evolution, the positions of the interfaces could suddenly be frozen, then the 
liquid region would subsequently relax exponentially towards a uniform gradient. 
(Compare the related diffusion problem of the temperature profile in a conducting 
rod with the ends held at different temperatures.) In this sense, the liquid region 
is constantly in a state of relaxation towards a uniform density gradient, for which 
equation [16] applies. 

Let us define A to be the size of the metastable liquid region; A = x\ —X2- Then 



substituting equations 15 and 16 in O and 14 gives us the differential equation 



d(t/r) w ' 2A 
with a constant 

7i ( p B - Pa\ 2 



" A Wr)-i - £ (17) 



which has units of time/ length 2 , and a dimensionless constant 7, which acts as an 
attractive coupling between the walls, given by 

^2.(e^) 2 ?i lfiB - Pc )(^ + ^-). (19) 

V & J Di \p B -p A PD-pcJ 

The formula for 7 is quite easy to understand. It is a ratio of quantities which drive 
the interfaces together (relating to diffusion in the metastable phase) to those which 
drive them apart (relating to the gaseous region). It is proportional to the ratio 
of the diffusion constants, and to the difference in densities (ps — Pc) which drives 
flux through the liquid region. This difference is made dimensionless by the factor 
following it, which is dominated by the interface of smallest height. These quantities 
are divided by the square of the relative supersaturation, which is responsible for 
driving the g-l wall away from the l-s wall. 

4.3 Solution 

In equation [17], 7 has a critical value of unity, above which the solution may be 
expressed parametrically as 

A . . 9 



t sin 9 exp — 

V 7 - 1 V 7 
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Ai ' (sin 9 + a/7 - 1 cos 9) exp = (20) 



for values of the parameter 9 in the range < 9 < n — arctan ^7 — 1. Here, 
Ao is the initial size of the metastable region. A critical value of 7 implies (from 
equation |l^) a critical value a c of the adjusted supersaturation. The condition 7 > 1 
corresponds to a' < cr c . We see that a graph of A versus \/i is an affine deformation 
of a logarithmic spiral, and the restricted domain of 9 gives a branch thereof in the 
first quadrant. Hence the metastable region has a finite lifetime, since A decays 
to zero at some positive value of t. We will refer to these solutions as 'diffusively 
bound', to distinguish from the tighter binding due to curvature energy in the Cahn- 
Hilliard model, which prevents a single g-s interface splitting into g-l and l-s parts 
[§]. Furthermore, there is no solution for A = 0. Hence, for a' < a c , the flux 
of condensation cannot, even momentarily, separate the g-l from the l-s wall, if 
they are initially together. The diffusively bound solutions for A as a function of 
t/r are plotted with dotted lines in figure |6]a for various values of 7, with A = 1. 
Notice that the gradients of the curves are infinite where they meet each axis. The 
infinite gradient at t = arises from the infinite flux of condensing material due 
to the delta singularity in the density at the base of the g-l interface. At A = 0, 
the metastable region decays infinitely fast since the density gradient (p B — pc)/A 
diverges. Of course, in a physical situation, the gradient of A(t) would be flattened 
in both instances, since the densities (pa, Pb) do not truly remain constant for a 
fast-moving interface. Once the interfacial separation collapses to zero, the model 
breaks down. In a real system, the interfaces would combine into a single, stable 
g-s interface, which would continue to advance. 



At the critical value 7 = 1, equation [H] has the parametric solution 



t/r = A o (0 + 1) exp0 
A = A o 0exp0 (21) 

and for 7 < 1 (i.e. above the critical supersaturation), the solution is 
yt/T = smn <p exp 



V 1 -7 V 1 -7 

A = (sinh0 + Jl — 7 cosh0) exp = (22) 

VI - 7 v Vl - 7 

with < cf) < 00 in each case. Graphs of A versus t/r, with A = 1, are plotted in 
figure |^a, with solid lines for 7 < 1, and a dashed line for 7 = 1. These solutions 
are not bound, i.e. A > for all positive t and A — ► 00 as t — ► 00. So, above the 
critical supersaturation, the flux of condensing material, although always dwindling, 
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Figure 6: (a) Size of the metastable region A (as given in equations |20], ^T] and |22|) 
against time t in units of the constant r, for various values of the coupling constant 
7. Diffusively bound solutions are marked with a dotted line, for 7 = 3, 7 = 2 
and 7 = 1.8. The marginal solution (7 = 1) is dashed, and unbound solutions are 
plotted with solid lines, for 7 = 0.8 and 0.1. In each case, the initial size of the 
metastable region, A is unity, (b) The same plots, produced by numerical solution, 



without the approximations, equations [15] and |16[ The same values of 7 are used as 
in (a). 



is sufficient to separate the g-l from the l—s interface, causing the metastable liquid 
phase to grow ad infinitum. This is even true in the extreme case where Ao = 0, 
when the solution becomes 



A= M + ^/l- 7 J yjt/r for 7 > 1. (23) 

Equation ^ is also the limit of equation^ as t —>■ 00, so the size of the metastable 
region at late times is independent of its initial value. 

To check the validity of equations [15| and [16|, we have solved numerically for 
the system described in section |4.1| (without the approximations made in section 
|4.2| ). The results plotted in figure |]b are for D2/D1 = 0.1. The variation of 7 was 
controlled by varying a' only. These results compare well with the approximate, 
closed-form solutions in figure |6]a. The same qualitative features (open/closed tra- 
jectories) appear, and the critical value of 7 is close to unity. The lifetimes of the 
bound solutions (which are very sensitive to the system parameters) agree to within 
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a few percent for systems not too close to criticality. For the top-most trajectory in 
figure ||b, a' / (pb~ Pa) ~ 50% so some deviation from linearity (equation |l^) is to be 
expected (see section |3lf ). We conclude that equations [L5| and [16] are quantitatively 
reasonable, and qualitatively very good, approximations. 

In summary, the asymptotic (late-time or long-distance) behaviour of a metastable 
region, for a system in which a double interface has formed early in the phase or- 
dering, is as follows. For an ambient supersaturation a' below the critical value 



0"c = (pb ~ Pa) 



\ 



27t(pb - Pc) 




(24) 



there is no asymptotic behaviour, since the metastable region collapses in a finite 
time, and subsequent evolution is via the ordinary g-s mode of interface propagation. 
Above the critical supersaturation, the metastable phase grows with time, according 
to equation |23[. Notice, with reference to figure |2|b, that (p B — p c ) is a measure of 
the metastability of the liquid phase. Hence at the triple point, when the three 
wells have a single common tangent, the critical supersaturation given by equation 
2~4] goes to zero. In other words, at the triple point, no supersaturation is required 
to stabilize the liquid phase, so the physics is modelled successfully in this respect. 

So far, we have concentrated on the size, A, of the metastable region, but the 
growth of the solid region is also of interest. For the double-interface mode of 



evolution, the size of the solid region is given, from equations [14| and [Rj, by 



x 2 {t) 



Pb - Pc 
,Pd - Pc, 



Do 



* dt' 
o A{F) 



Substituting the asymptotic expression for A (equation 
using equation 



into this formula, and 



18 for r, gives the result 



J' 2 



2v^F 



Do 



Pb - Pc\ f Pb - Pa 



(1 + V 1 - 7) \fD~i \Pd ~ Pc, 



a 1 



Notice that this is inversely proportional to the (adjusted) supersaturation. It is 
interesting to compare this with the size of a solid region produced by normal inter- 
face motion (i.e. by a single g-s interface). To compare like with like, we should use 
the linearized velocity coefficient (equation [IT], with (p B — Pa) replaced by (p s — p g ) 
for the interface height), together with equation |8|. Denoting the size of the solid 
region produced by single interface motion by x S i ng ., we find 



X s 
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Let ( be the ratio of the size of the solid region produced by the double-interface 
mode of growth, x 2 , to that produced by normal growth, x sing .. Then 

71 (PB ~ PA\ ( Ps - Pg\ D 2 {p B ~ Pc) 



C = 71 ( PB~ PA \ ( Ps- Pg \ 



V 1 ~ 7) V cr' / V a J Di (p D - Pc) ' 
There is a distinct similarity between this formula, and the expression for 7 (equation 



19|) . If we approximate (pp — Pa) + (pb — Pc) by (p s — p g ), which we see, from figure 

c 



0b, is usually a good approximation, then we find 

7 a 



2(1 + vT=t) a ' 

Since < 7 < 1 for the split-interface mode of growth, and also a' < cr, it follows 
that 

C < ~ (25) 

and that, well above the critical supersaturation, £ w 7/4. 

Let us recapitulate the properties of the parameter 7. It appeared in equation 



17| as an attractive coupling between the g-l and l-s interfaces. It is the ratio of 
properties of the liquid phase, tending to attract the interfaces, to properties of the 
gaseous phase, tending to separate them. It characterizes the classes of solutions of 
equation O, being greater than unity for diffusively bound solutions, and less than 



unity for unbound solutions. Finally, we see in equation |25] that it gives an upper 
bound (and an order of magnitude) for the ratio of growth rates of the solid phase 
in the split and normal modes of growth. Since 7 < 1 when the split mode occurs, 
equation |25] shows that this mode of evolution suppresses the growth of the solid. 

Note that, in figure §b, the metastability of the middle well leads to the inequality 
Pb > Pc- If this liquid well were stable (i.e. below the double tangent to the outer 
two wells), the inequality would be violated. Hence it is reasonable to define the 
metatability, m, of the middle well to be m = ps — pc- With reference to equation 
|23|, we see that the boundary between regions of the (m, <r')-plane, for which the g-l 
and l-s interfaces are diffusively bound or unbound, is of the form shown in figure 0. 
A naive expectation would be for a boundary coincident with the vertical axis, but 
we see that this is not the case. (For negative a', for which the liquid phase must 
dissolve into the gas, we find solutions of equation |17| are again given by equations 



20| , pi| and but with different ranges of the parameters 9, and (ft, leading to closed 



trajectories whenever 7 > 0.) 

5 Discussion 

Some non-equilibrium effects in colloid-polymer mixtures have been reviewed in 
section |2|. We have remarked on the presence of a well-defined non-equilibrium 
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Figure 7: The regions of metastability m = p B — p c and adjusted supersaturation 
a' for which interfaces are diffusively bound or unbound. 



boundary in the phase diagram, close to the theoretically calculated position of 
the 'hidden' gas-liquid metastable binodal. This metastable binodal is central to 
the theoretical model developed in section |], where the position of its low-density 
branch was denoted by pa- A system homogenized by shear melting to an ambient 
density of exactly pa would, in the notation of section |] have an vanishing 'adjusted 
supersaturation', a' = 0. In the model, this value has the special significance that, 
for a' > 0, a region of metastable liquid has an initial period of growth, before it 
collapses. At lower densities, any liquid region that is nucleated will immediately 
shrink. In experimental samples with a colloidal density above the non-equilibrium 
line in the phase diagram, the growth of crystals is observed to be suppressed. This 
happens in the model for a' > a c . 

In practice the lines o' = (the metastable binodal) and a' = a c may be ex- 
perimentally indistinguishable, for two reasons. Firstly, no very accurate theoretical 
prediction exists for the position of the metastable binodal in the phase diagram. 
Attempting to make quantitative comparisons of time scales with experiment, using 
the best available theories for the values of pa, Pb, Pc an d Pd, results in uncertain- 
ties of several hundred percent in the value of 7. Secondly, the line a' = a c may be 
very close to the metastable binodal if a c is small. (From equation [24], this would 
occur if e.g. the middle well in the potential were only slightly metastable, or if the 
diffusion constant were much lower in the liquid than in the gas.) If a c is small 
then, writing 7 = {a c /a') 2 for p B > p c , we see that 7 will decrease rapidly with 
increasing colloidal density, and so, from equation |25, suppression of crystal growth 
will be pronounced, at densities not far from the metastable binodal. 
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It seems, then, that the model presented here provides a plausible theoretical 
basis for the previous conjecture that the onset of nonequilibrium behaviour in 
certain colloid-polymer mixtures is associated with the presence of a 'hidden' gas- 
liquid binodal. It is, however, a greatly simplified and idealized picture, and we 
should consider the ways in which it deviates from reality, and the implications for 
the resultant interfacial dynamics. 

The model is one-dimensional, and therefore ignores surface tension. This is jus- 
tifiable since, once a region has grown considerably larger than the critical nucleus 
size, surface tension has a negligible effect on interfacial motion during the growth 
stage [j!4| . It is at these intermediate-to- late times that our model describes the 



system. Dimensionality is also relevant to the time dependence of the long-range 
diffusion. This will have a quantitative effect on the predicted values of the critical 
supersaturation and the degree of suppression of crystal growth etc., but we may 
conjecture that the qualitative features of the model's behaviour will extend to three 
dimensions. The fact that our one-dimensional model does not explicitly address 
intrinsically higher- dimensional geometric effects, such as the Mullins-Sekerka insta- 
bility [IS|, may not be important. The model gives us the general rule of thumb that, 
above a certain critical supersaturation, the crystalline regions (whatever shapes 
they may be), which would normally grow in a two- well system, are replaced, in 
the presence of a third well, by metastable liquid. Subsequent to this growth stage, 
the metastable liquid is slowly transmuted into the 'correct' equilibrium phase i.e. 
crystal. 

Probably a more drastic simplification is the semi-infinite extent of the 'ambient' 
gaseous region in the model. In reality, there is more than one nucleus of the solid 
phase. There will be some typical inter-nucleus spacing L nuc in the system. The 
concentration profile in the gaseous region has a characteristic length scale y/(Dit) 
so, when this is of order L nuc ., the nuclei begin to influence each other. After this 
time, the effective asymptotic supersaturation begins to fall, as the regions of deple- 
tion of the gas phase around the nuclei begin to overlap. This characteristic time to 
deplete the supersaturation of the gaseous phase may be denoted tde P .(~ L^ uc /D^. 
At this point, the metastable region around each nucleus is of size A(td ep .), which 
is given in equation The remaining lifetime of the metastable regions may be 
calculated from equations |2D], as the time for a region initially of this size to col- 
lapse, in a system with adjusted supersaturation of zero. The total lifetime of the 
metastable phase is thus of order idep./7- (Compare this with the naive calculation 
for a semi-infinite gaseous region, which gives a suppression of size of the solid region 
by a factor of order 7, and hence a time-scale factored by I/7 2 , rather than 1/7.) 
We note that, if a c is small, then the lifetime of the metastable phase, (a'/a c ) 2 td ep . 
grows rapidly with the ambient density. Once the metastable regions have collapsed, 



20 



and g-s interfaces form, the densities will diffusively readjust from pa and po to p g 
and p s , in a relatively short time. 

This whole discussion assumes that all interfaces in the system are undergoing 
the 'split' mode of propagation (i.e. with g-l and l-s parts not bound together by 
curvature, as discussed in section || and reference ||). Any "normal" (g-s) interfaces 
initiated in the system during nucleation will lead to the formation of large crys- 
tals on normal diffusive time-scales. Such crystals are not observed experimentally 
above the non-equilibrium boundary. Therefore, our model is a good candidate for 
the physics of the non-equilibrium boundary if, for some reason, only split inter- 
faces are generated during nucleation above this boundary. Such a scenario is not 
unreasonable. We have seen that diffusively bound interfaces are unbound by a 
sufficiently large supersaturation. It seems likely || (and has in fact been observed 
in a preliminary numerical study ||) that a g-s interface (which is stabilized, or 
"bound" by a contribution to the chemical potential of the form — V 2 p) may be 
"split" or "unbound" in a manner analogous to diffusive unbinding and, further- 
more, that the critical supersaturation to cause this is close in magnitude to that 
calculated here. This conjecture is based on the fact that, in a model that includes 
an extra V 2 p term in the chemical potential to describe curvature effects (such as 
the Cahn-Hilliard model), the curvature term exactly balances the diffusive term in 
an equilibrium interface and hence quantities calculated from it will, in the main, 
be of the same order of magnitude as those calculated from the diffusive term only. 



6 Related systems 

The model we have investigated in this paper was originally suggested by experi- 
mental observations in mixtures of spherical colloids and non-adsorbing polymers of 
a substantially smaller size. The same model may, however, be applicable to other 
experimental systems. 

First of all, it has been pointed out recently by one of us that the 'hidden binodal' 
is probably significant for understanding the crystallization of globular proteins |16| . 



The kinetic predictions of this paper may therefore also be relevant in that context. 

Furthermore, our model is probably directly relevant to mixtures of rod-like col- 
loids and non- adsorbing polymers [17|. The two possible phase diagram topologies 



for this system are again those given in figures ||a and ||b, but with the labels gas, 
liquid and solid replaced by 1%, I2 and iV (standing for isotropic phases 1 and 2, 
and nematic phase). Once more, three-phase coexistence disappears when the size 
of the polymer (relative now to the rod length) decreases below a critical value. In 
this case (compare figure ||b), addition of sufficient polymer to a suspension of rods 
leads to slow phase separation into an isotropic and a denser, coexisting nematic 
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phase, the latter being distinguished by strong birefringence. Further addition of 
polymers, however, brings about a different kind of behaviour — a weakly birefrin- 
gent, 'expanded' phase, which contains most of the rod-like particles, separates out 
quickly |18| . It has been suggested [18| that the suppression of isotropic-nematic 



phase separation is due to the presence of a 'hidden' isotropic-isotropic binodal. 

Let us return to spherical colloid and polymer mixtures. In the preceding sec- 
tions, we have concentrated on polymer-colloid size ratios £ sufficiently small for no 
liquid phase to appear in the equilibrium phase diagram (figure |l|b). Now consider 
the case where the polymer is large enough to give rise to a thermodynamically sta- 
ble gas-liquid binodal, as appears in figure pi. Note that, in the 'gas+solid' region 
of figure |l|a, the form of f(p) is as sketched in figure §b, and hence the model set 
up in this paper is again relevant. An expanded version of the phase diagram near 
the triple line, now showing the metastable portions of the gas-liquid binodal, is 
sketched in figure [| As the triple line is approached from above (p — > p P r ), where 
the liquid well is only just metastable, the critical supersaturation vanishes (cr c — > 0) 
because (p B — p c ) — > 0; see figure ^]b and equation ^4} At higher values of p p , we 
expect another regime where a c — > 0, this time due to the vanishing of the diffusion 
coefficient in the liquid phase, D 2 (see equation |24[). As p p increases, the density of 
the metastable liquid phase, ps also increases; eventually ps — > Pgiass, the density 
at which the system vitrifies, with D 2 — > 0. 

We can thus speculate on the form of the boundary for diffusive unbinding, 
p = pa(p p )+ct c (pp), in the vicinity of the gas branch of the metastable binodal, figure 
H The prediction is that, for moderate colloid densities, there should be a region 
of suppressed crystallization immediately above the triple coexistence line (due to 
the very marginal metastability of the liquid minimum in the free energy density), 
followed, for higher p p , by normal crystallization behaviour, and ending up with 
crystallization suppression again at the highest values of p p (due to the vitrification 
of the metastable liquid phase). Experimentally, the former non-equilibrium region 
has not been observed fllTf . A search for this phenomenon is under way in our 



laboratory. However, the region of normal crystallization has been reported, as has 
a non-equilibrium boundary at higher p p where it has been suggested [|19j that 
vitrification of a dense phase does play a crucial role. 

We speculate that the theoretical results contained in this paper should have 
some relevance to a number of other complex fluid systems in which metastability is 
known to play a key role in phase transformations, including the crystallization of 
the 'monotropic liquid crystal', poly-n-nonyl-4 4'-biphenyl-2-cholroethane, via an in- 
termediate, metastable nematic phase |HJ and the crystallization of poly(phenylene 
ether) in cyclohexanol, where deep quenches produces first a fluid-fluid phase sepa- 
ration |21| . We acknowledge, however, that the limitations of our model discussed 
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in section [5], together with the likelihood that latent heat may not be negligible in at 
least some instances, necessitates further research before reaching firm conclusions 
on its applicability. 




P 

Figure 8: The phase diagram in figure |l|a redrawn schematically to emphasize 
the portion above the triple coexistence line at p} p r . Here, within the equilibrium 
gas+solid coexistence region, the metastable gas-liquid binodal is the dotted lines 
p = Pa(h p ) and p = p#(/ip). The bold curves indicate the likely positions of parts of 
the diffusive unbinding boundary p = pa + Cc(p P )- Near the triple coexistence line, 
the boundary meets the gas branch of the metastable binodal because the liquid 
well in the free energy density is only marginally metastable. At pf ass , the liquid 
branch of the metastable binodal reaches the vitrification density, p g i aS s- in this 
vicinity, the critical supersaturation curve is again close to the gas branch of the 
metastable binodal, now because of the vanishing diffusion constant in the liquid 
phase. NC indicates regions where we expect disruption of crystallization, while 
normal gas-crystal coexistence is expected in the hatched area. 



7 Conclusions 

We have studied a simple model for diffusion-limited kinetics of phase ordering 
in a system whose free energy density has a metastable third well, at a density 
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intermediate to the two equilibrium phases. In such systems, we have found a 
mechanism whereby a region of the metastable phase may grow ad infinitum at the 
expense of the equilibrium dense phase, if the mean density of the system is above a 
critical threshold. This behaviour appears to be consistent with the non-equilibrium 
ordering dynamics and suppression of crystal growth observed experimentally in 
colloid-polymer mixtures. In the experimental system, as in the theoretical model, 
the onset of anomalous behaviour occurs at a well defined density. 

The important lesson of this study (and the companion paper ||) is that any 
pair of concentrations, which can be linked by a double-tangent on the graph of free 
energy density versus concentration, may give rise to an interface in the evolving 
system. Although only the globally stable binodal densities will coexist in the equi- 
librated system, local and transient coexistence can occur between any binodal pairs 
of densities in a system which has not yet discovered its global equilibrium state. 
Hence, metastable phases cannot be overlooked when modelling phase ordering. In- 
deed, their importance has long been accepted on an empirical level, particularly 
in a metallurgical context [E2L [TBI]. Ostwald's 'rule of stages', for example 



asserts that the transformation from one stable phase to another proceeds via all 
metastable intermediates in turn. We have, in the present paper, a rationale for the 
consideration of such phases in a soft condensed matter setting. 
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A Evaluation of Integrals 

Dividing equation | by equation [7], and using the appropriate functional forms for 
Xi(t) and s(t), and the substitution t' = u 2 t to eliminate dimensionality from the 
integrands, we find the relation between the velocity coefficient and the relative 
supersaturation <j/(pb — Pa) = h/h where 

f 1 du a 2 (l + u) 

Jo exp ~ 4 (l- M ) 
r 1 du a 2 (l + u) 

Jo (1 - u)VT^? 6XP ~ 4(1 - u) ' 



h = 
and I2 = 
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Making a further change of variables v 2 



(1 + «)/(! -it) yields 



and I\ 



f 




dv. 



We see that I\ and I2 are related by a derivative: 



da 



[he-* 



2 




Hence I\ is given by the indefinite integral 



h 




I 



e 




whose constant of integration is set by noting that I\ vanishes as a — >• 00. This 
integral is soluble by parts, using 



from which equation |9| follows. 
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